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Abstract 

The mechanism of the flux-line-lattice (FLL) melting in anisotropic high- 
ly superconductors in B || c is clarified by Monte Carlo simulations of the 
3D frustrated XY model. The percentage of entangled flux lines abruptly 
changes at the melting temperature T m , while no sharp change can be found 
in the number and size distribution of vortex loops around T m . Therefore, 
the origin of this melting transition is the entanglement of flux lines. Scaling 
behaviors of physical quantities are consistent with the above mechanism of 
the FLL melting. The Lindemann number is also evaluated without any 

phenomenological arguments. 
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Nature of the mixed phase in type-II superconductors has been studied for many years, 
and much attention has been paid to this field since the discovery of high-T c superconductors 
(HTSC) because of short correlation lengths and large anisotropy. In HTSC in a magnetic 
field along the c axis, the flux-line lattice (FLL) melts at much lower temperatures or in 
much weaker fields than those predicted by the Abrikosov mean-field theory, where the 
superconducting phase transition is of second order regardless of details of models. The FLL 
melting in HTSC was first theoretically analyzed by Nelson and Seung H on the basis of 
the mapping to a two-dimensional Boson model, and they pointed out that the Lindemann 
criterion might be valid in this FLL melting because of large fluctuations in HTSC owing 
to large anisotropy and high transition temperatures. Similar theoretical analysis was also 
made by Houghton et al. [Q] independently. Earlier experiments of "the FLL melting in 
HTSC" as reviewed in Ref. 1 were found to be explained better by the vortex-glass transition 
0] rather than by the FLL melting transition. In clean systems, the "true" FLL melting 
was confirmed afterwards by experiments |5|-|T0f and computer simulations [p"l|-|T8|j, and the 



first-order FLL melting transition in a magnetic field has now been established in extremely 
type-II superconductors such as HTSC. 

On the other hand, the mechanism of the FLL melting has not yet been well understood. 
Nelson argued Pl,|l9| that the entanglement of flux lines is related to this transition, but 
this picture has only been confirmed numerically |l7j within the two-dimensional Boson 



model. Thermal excitations of vortex loops [^,^1 are n °t included in this picture and some 
authors claimed plt|22| that the FLL melting is characterized by the breakdown of flux-line 
description by the proliferation of large vortex loops, but the latter picture began to be 
modified recently f23|| . In order to clarify the mechanism of the FLL melting numerically, 



the first-order phase transition should be identified by measuring thermodynamic quantities, 
and the behavior of flux lines induced by an external magnetic field and thermally excited 
vortex loops should be observed microscopically in the vicinity of the melting temperature. 

In this article, the three-dimensional anisotropic, frustrated XY model is analyzed with 
the Monte Carlo method from the above point of view. Our main results are as follows: First, 



the mechanism of the first-order FLL melting transition is exclusively the entanglement of 
flux lines. The percentage of entangled flux lines sharply changes at the melting temperature 
T m , while the number and the size distribution of loop excitations has a smooth temperature 
dependence around T m . Second, scaling properties around the melting temperature are 
clarified. As a consequence of the entanglement mechanism of the FLL melting, T m is scaled 
by the inverse of the system size along the c axis. Third, the Lindemann number takes 
nearly a constant value cl ~ 0.30 regardless of the anisotropy constant, and thus the use 
of the Lindemann criterion is justified for the determination of the melting line in a phase 
diagram. 

As the model of the anisotropic, extremely type-II HTSC in a magnetic field along the c 



axis, we consider the three-dimensional anisotropic, frustrated XY model p^Jlql described 
by the following Hamiltonian, 

H = —J Yj cos (<A - Vj ~ A ij) 

i,j£ab piano 
J 



■ 2 COs((Pi~<Pj) , (1) 



i,j\\c axis 

4* = X / J A«.dr« (2) 

00 Ji 

Here (pi denotes the phase of the superconducting order parameter, O stands for the flux 
quantum, and the anisotropy is represented by the parameter 1/T 2 . Note that this model 
neglects fluctuations of the gauge field and the amplitude of the superconducting order 
parameter. The periodic boundary condition (PBC) is applied on the phase variable ifii in 
all the directions in order to refrain from finite-size effects from free boundaries. We mainly 
investigate the case with the averaged number of fluxes per plaquette / = 1/25,T = 2 
and 5. In these parameters, effects of the introduction of a square lattice in the ab plane is 



negligible [15]. Each phase variable takes a value — n < <pi < n, and the summation of the 



phase difference around a plaquette is given by 



5^ (<Pi - <pj - Aij ) = 27r(n - /) , n = 1, 0, -1 . (3) 

i,jen 



When the integer n takes 1 or —1, the plaquette is defined to have a vortex or an antivortex, 
respectively. The nearest-neighbor vortices are connected with one another to form vortex 
lines, which do not have end points inside the system. When a vortex line returns to itself 
inside the system, it is called as a vortex loop, and such a loop does not exist in the ground 
state. When a vortex line runs from one boundary to another along the direction of the 
external field, it is called as a flux line. The flux lines are straight in the ground state, 
and they begin to fluctuate at finite temperatures. A flux line is defined as entangled if it 
does not terminate at the same transverse position in the top and bottom boundaries in the 
PBC. The flux lines which wind with each other and return to the initial transverse positions 
inside the system are not included in this definition. Such excitations are negligible in the 
vicinity of the melting temperature, because they have higher energies. 

The helicity modulus along the c axis |T5| , |25[| is observed for the determination of the 
melting temperature. This quantity is proportional to the superfluid density, and therefore 
nonvanishing only in the superconducting phase. The numbers of vortex loops A^ oop and 
entangled flux lines N ent are counted. The distribution of sizes is also measured for Josephson 
loops, which are dominant below and slightly above the melting temperature. The transverse 
distance w of a flux line between the top and bottom ab planes is measured, and its averaged 
value over all the flux lines is denoted by Ldiff • The fluctuation of a flux line is measured by 
the deviation u from the projection of its mass center in each ab plane, and averaged over 
all the flux lines and the ab planes. The Lindemann number cl is defined by 

c L = lim (uY 2 /a Q , (4) 

1 > 1 m U 

where ao stands for the lattice constant of the triangular FLL, ao = (2/v / 3) 1 ^ 2 // 1 ^ 2 - 

Monte Carlo simulations are performed on the basis of the Metropolis algorithm. Most 
results reported in this article are for systems with L x = L y = 50 and L c = 80. Then, the 
number of total flux lines is iVfl ux = 100 in the present simulations with / = 1/25. In order 
to check the size dependence, we also simulate systems with (L X ,L C ) = (50,20), (50,40), 
(50,54), (50,160), (25,80) and (100,20). In addition, in order to check the flux-density 
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dependence, the (L X ,L C ) = (50,80) and (100,80) systems are calculated for / = 1/50 and 
1/100, respectively. Simulations are started from temperatures more than ten times higher 
than T m and the system is gradually cooled down. Typical Monte Carlo steps (MCS) are 
1.0 x 10 5 and 1.5 x 10 5 for equilibration (E-MCS) and measurement (M-MCS) at each 
temperature, respectively. Since the correlation time becomes longer in the vicinity of the 
melting point, E-MCS and M-MCS are taken as 3.5 x 10 5 and 8.0 x 10 5 , respectively, and the 
cooling rate is reduced to AT = 1.0 x 10 _3 J//cb- Moreover, for the precise determination 
of T m and cl, supercooling behavior is reduced very carefully by using up to 1.0 x 10 7 MCS 
at each temperature. The helicity modulus is measured at each MCS, and the numbers of 
vortex loops and entangled flux lines are measured once per 100 MCS. 

The temperature dependence of the helicity modulus along the c axis, T c , is displayed 
in Fig. |l] for T = 2 and 5 with / = 1/25. This quantity sharply drops from a finite value 
to zero at the melting temperature, T m ~ 0.810 J/k-& for Y = 2 and T m ~ 0.3445 J/k-Q for 
T = 5, which indicates the thermodynamic first-order phase transition [HJ. The temperature 
dependence of the ratio of entangled flux lines to total flux lines, iV ent /iVfl ux , is also displayed 
in Fig. [TJ. It shows a sharp jump at T m for each of the anisotropy. Similar behavior is also 
observed for V = 5 with / = 1/50 and 1/100. The number of vortex loops per flux line per ab 
plane, iV"i oop /(./Vfl ux X c ), is shown in Fig. § for T = 2 and 5 with f = 1/25. The temperature 
dependence of this quantity is not as drastic as that of the ratio of entangled flux lines. The 
size distribution of Josephson loops is also measured for T = 5 with / = 1/25, 1/50 and 
1/100, and no drastic change is observed in this distribution around T m . The numbers of 
vortex loops are not of the same order for the different anisotropy constants at the melting 
temperatures. That is, the number of vortex loops at T m for T = 5 corresponds to that 
of T « 0.6J//cb for T = 2, and this temperature is much lower than the melting point for 
T = 2, T m ~ 0.810 J/kB- These facts clearly indicate that the origin of the FLL melting is 
the entanglement of flux lines, at least up to / = 1/100. 

Then, we show the results for the finite-size-scaling behavior. The averaged end-to-end 
transverse distance of flux lines, L^m, is normalized by the lattice constant of FLL, ao, and 
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plotted versus temperature for L c = 40, 80 and 160 in Fig. ||. The size dependence of this 
quantity can be described by the random-walk-type scaling, 

L d is(L c ) ~ const, x L\ /2 , (5) 

for a certain temperature range above T m , as displayed in the inset of Fig. |3|. Therefore, it is 
clear that the vortices form flux lines even above T m . The inset of Fig. |3| shows that the data 
for L c = 160 and 80 deviate from the scaling at T « OAJ/ks and 0.5 J/ ks, respectively. 
These two temperatures correspond to the same transverse distance Ldiff ~ 2.7ao, as can 
be read from the two curves for L c = 160 and 80 in the main body of Fig. |3|. This fact 
suggests that the random-walk behavior is restricted within a transverse diffusion distance 
approximately 2.1 independently of temperature even in bulk systems. Beyond this length 
scale, reconnection between flux lines occurs frequently, and the random-walk property is 
suppressed. Since the transverse distance between the bottom and top ab planes in each 
entangled flux line cannot be smaller than a , the total number of entangled flux lines 
abruptly decreases at T m when temperature is gradually reduced. As a consequence of the 
scaling given in Eq. (|5|), the temperature characterized by L diS ps a , namely T m , depends 
on L c , as will be discussed below. 

Simulations for systems with (L X ,L C ) = (100,20) and (25,80) are also performed, and 
the results coincide with those with (L X ,L C ) = (50,20) and (50,80), respectively. This is 
quite natural because the melting transition is characterized by the entanglement of flux 
lines along the c axis, and the leading term of size dependence is only related to L c . Finite- 
size effects in the ab plane are indirect on thermodynamic quantities. As displayed in Fig. |], 
our data exhibit the following finite-size scaling, 

5T m (L c ) = T m (L c ) - T m (L c = oo) ~ const, x L' 1 , (6) 

with T m (L c = oo) = 0.3354 ± 0.0007 for T = 5. This size dependence means that the FLL 
melting has a one-dimensional character, because the scaling form of the transition point in 
first-order phase transitions [215 is generally given by ST(L) ~ const, x L~ D with the spatial 



dimension D. This one-dimensional character is consistent with the entanglement picture. 
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Finally, we turn to see the temperature dependence of the fluctuation of flux lines and 
evaluate the Lindemann number. The quantity (M 2 ) 1//2 /a shows a sharp jump at T m for 
T = 2 and 5 with / = 1/25 (Fig. ^]). From the definition of the Lindemann number in 
Eq. (f|), we have cl ~ 0.30 for both the anisotropy constants. This result suggests that the 
Lindemann number does not depend on details of models, as assumed in previous studies 
|3|j27[|. The Lindemann number was evaluated as cl ~ 0.18 by part of the present authors 
|15| by fitting the simulated melting line with a formula derived by Blatter et al. [27 based 
on the London theory. We believe that the present direct evaluation of cl is more reliable. 

In conclusion, the three-dimensional anisotropic, frustrated XY model has been analyzed 
with Monte Carlo simulations. The melting temperature T m has been estimated as the 
point at which the helicity modulus along the c axis vanishes. The percentage of entangled 
flux lines shows a sharp jump at T m , while the number and size distribution of vortex 
loops do not show such drastic change at T m . This fact clearly indicates that the origin 
of the FLL melting in a magnetic field along the c axis is the entanglement of flux lines. 
The consistency of this picture with the size dependence of various quantities has been 
confirmed. Especially, the melting temperature is scaled by the system size along the c 
axis as T m (L c ) — T m (L c = oo) oc L^ 1 . The averaged deviation of flux lines from their 
mass centers also shows a sharp as a consequence of the entanglement of flux 

lines. The Lindemann number takes a constant value cl ~ 0.30 regardless of the anisotropy. 
This numerical result justifies the use of the Lindemann criterion for characterizing the FLL 
melting in HTSC. 

Numerical calculations were performed on the Numerical Materials Simulator (NEC SX- 
4) at National Research Institute for Metals, Japan. 



7 



REFERENCES 

[1] G. Blatter et al, Rev. Mod. Phys. 66, 1125 (1994). 

[2] D. R. Nelson and H. S. Seung, Phys. Rev. B 39, 9153 (1989). 

[3] A. Houghton, R. A. Pelcovits, and A. Sudb0, Phys. Rev. B 40, 6763 (1989). 

[4] D. S. Fisher, M. P. A. Fisher, and D. A. Huse, Phys. Rev. B 43, 130 (1991). 

[5] H. Safar et al, Phys. Rev. Lett. 69, 824 (1992). 

[6] W. K. Kwok et al, Phys. Rev. Lett. 72, 1092 (1994). 

[7] U. Welp et al, Phys. Rev. Lett. 76, 4809 (1996). 

[8] J. A. Fendrich et al, Phys. Rev. Lett. 77, 2073 (1996). 

[9] E. Zeldov et al, Nature (London) 375, 373 (1995). 
[10] A. Schilling et al, Nature (London) 382, 791 (1996); Phys. Rev. Lett. 78, 4833 (1997). 
[11] R. E. Hetzel, A. Sudb0, and D. A. Huse, Phys. Rev. Lett. 69, 518 (1992). 
[12] R. Sasik and D. Stroud, Phys. Rev. Lett. 75, 2582 (1995). 
[13] S. Ryu and D. Stroud, Phys. Rev. B 54, 1320 (1996). 
[14] S. Ryu and D. Stroud, Phys. Rev. Lett. 78, 4629 (1997). 

[15] X. Hu, S. Miyashita, and M. Tachiki, Physica C 282-287, 2057 (1997); Phys. Rev. Lett. 
79, 3498 (1997); Phys. Rev. B 58, 3438 (1998). 

[16] J. Hu and A. H. MacDonald, Phys. Rev. B 56, 2788 (1997). 

[17] H. Nordborg and G. Blatter, Phys. Rev. Lett. 79, 1925 (1997); Phys. Rev. B 58, 14556 
(1998). 

[18] A. E. Koshelev, Phys. Rev. B 56, 11201 (1997). 



8 



[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 



D. R. Nelson, Phys. Rev. Lett. 60, 1973 (1988). 
Z. Tesanovic, Phys. Rev. B 51, 16204 (1995). 

A. K. Nguyen, A. Sudb0, and R. E. Hetzel, Phys. Rev. Lett. 77, 1592 (1996). 
A. K. Nguyen and A. Sudb0, Phys. Rev. B 57, 3123 (1998); 58, 2802 (1998). 
Z. Tesanovic, |cond-mat / 980 YdUE . 



Y. -H. Li and S. Teitel, Phys. Rev. Lett. 66, 3301 (1991); Phys. Rev. B 47, 359 (1993). 

M. E. Fisher, M. N. Barber, and D. Jasnow, Phys. Rev. A 8, 1111 (1973). 

For example, D. P. Landau, in Finite Size Scaling and Numerical Simulation of Statis- 
tical Systems, ed. V. Privman (World Scientific, Singapore, 1990), p. 223. 



[27] G. Blatter, V. Geshkenbein, A. Larkin, and H. Nordborg, Phys. Rev. B 54, 72 (1996). 



9 



FIGURES 

FIG. 1. Helicity modulus along the c axis (squares) and the ratio of entangled flux lines to 
total flux lines (circles) versus temperature for T = 2 and 5 with / = 1/25. 

FIG. 2. Normalized number of vortex loops (diamonds) versus temperature for T = 2 and 5 
with / = 1/25. The ratio of entangled flux lines (circles) is also plotted for comparison. 

FIG. 3. Normalized end-to-end transverse distance versus temperature for L c = 40 (triangles), 
80 (squares) and 160 (circles) for T = 5 with / = 1/25. Scaling plot of the same data according to 
Eq. (||) is shown in the inset with the same symbols. 

FIG. 4. L c dependence of the melting temperature for T = 5 with / = 1/25. T m (L c = oo) is 
estimated by the least-squares fitting of the data with L c = 40, 54, 80 and 160. 

FIG. 5. Fluctuation of flux lines (triangles) versus temperature for T = 2 and 5 with / = 1/25. 
The helicity modulus along the c axis (squares) is also plotted for comparison. 
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Fig.1: Y.Nonomura, X.hu, and M.Tachiki 
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Fig. 2: Y.Nonomura, X.Hu, and M.Tachiki 
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Fig. 3: Y.Nonomura, X.Hu, and M.Tachiki 




Fig. 4: Y.Nonomura, X.Hu, and M.Tachiki 




Fig. 5: Y.Nonomura, X.Hu, and M.Tachiki 




